3.10.23 \(\int x^2 \sqrt {c x^2} (a+b x)^n \, dx\) [923]

3.10.23.1 Optimal result
3.10.23.2 Mathematica [A] (verified)
3.10.23.3 Rubi [A] (verified)
3.10.23.4 Maple [A] (verified)
3.10.23.5 Fricas [A] (verification not implemented)
3.10.23.6 Sympy [F]
3.10.23.7 Maxima [A] (verification not implemented)
3.10.23.8 Giac [B] (verification not implemented)
3.10.23.9 Mupad [B] (verification not implemented)

3.10.23.1 Optimal result

Integrand size = 20, antiderivative size = 131 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=-\frac {a^3 \sqrt {c x^2} (a+b x)^{1+n}}{b^4 (1+n) x}+\frac {3 a^2 \sqrt {c x^2} (a+b x)^{2+n}}{b^4 (2+n) x}-\frac {3 a \sqrt {c x^2} (a+b x)^{3+n}}{b^4 (3+n) x}+\frac {\sqrt {c x^2} (a+b x)^{4+n}}{b^4 (4+n) x} \]

output
-a^3*(b*x+a)^(1+n)*(c*x^2)^(1/2)/b^4/(1+n)/x+3*a^2*(b*x+a)^(2+n)*(c*x^2)^( 
1/2)/b^4/(2+n)/x-3*a*(b*x+a)^(3+n)*(c*x^2)^(1/2)/b^4/(3+n)/x+(b*x+a)^(4+n) 
*(c*x^2)^(1/2)/b^4/(4+n)/x
 
3.10.23.2 Mathematica [A] (verified)

Time = 0.05 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.74 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=\frac {c x (a+b x)^{1+n} \left (-6 a^3+6 a^2 b (1+n) x-3 a b^2 \left (2+3 n+n^2\right ) x^2+b^3 \left (6+11 n+6 n^2+n^3\right ) x^3\right )}{b^4 (1+n) (2+n) (3+n) (4+n) \sqrt {c x^2}} \]

input
Integrate[x^2*Sqrt[c*x^2]*(a + b*x)^n,x]
 
output
(c*x*(a + b*x)^(1 + n)*(-6*a^3 + 6*a^2*b*(1 + n)*x - 3*a*b^2*(2 + 3*n + n^ 
2)*x^2 + b^3*(6 + 11*n + 6*n^2 + n^3)*x^3))/(b^4*(1 + n)*(2 + n)*(3 + n)*( 
4 + n)*Sqrt[c*x^2])
 
3.10.23.3 Rubi [A] (verified)

Time = 0.21 (sec) , antiderivative size = 96, normalized size of antiderivative = 0.73, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.150, Rules used = {30, 53, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^2 \sqrt {c x^2} (a+b x)^n \, dx\)

\(\Big \downarrow \) 30

\(\displaystyle \frac {\sqrt {c x^2} \int x^3 (a+b x)^ndx}{x}\)

\(\Big \downarrow \) 53

\(\displaystyle \frac {\sqrt {c x^2} \int \left (-\frac {a^3 (a+b x)^n}{b^3}+\frac {3 a^2 (a+b x)^{n+1}}{b^3}-\frac {3 a (a+b x)^{n+2}}{b^3}+\frac {(a+b x)^{n+3}}{b^3}\right )dx}{x}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\sqrt {c x^2} \left (-\frac {a^3 (a+b x)^{n+1}}{b^4 (n+1)}+\frac {3 a^2 (a+b x)^{n+2}}{b^4 (n+2)}-\frac {3 a (a+b x)^{n+3}}{b^4 (n+3)}+\frac {(a+b x)^{n+4}}{b^4 (n+4)}\right )}{x}\)

input
Int[x^2*Sqrt[c*x^2]*(a + b*x)^n,x]
 
output
(Sqrt[c*x^2]*(-((a^3*(a + b*x)^(1 + n))/(b^4*(1 + n))) + (3*a^2*(a + b*x)^ 
(2 + n))/(b^4*(2 + n)) - (3*a*(a + b*x)^(3 + n))/(b^4*(3 + n)) + (a + b*x) 
^(4 + n)/(b^4*(4 + n))))/x
 

3.10.23.3.1 Defintions of rubi rules used

rule 30
Int[(u_.)*((a_.)*(x_))^(m_.)*((b_.)*(x_)^(i_.))^(p_), x_Symbol] :> Simp[b^I 
ntPart[p]*((b*x^i)^FracPart[p]/(a^(i*IntPart[p])*(a*x)^(i*FracPart[p]))) 
Int[u*(a*x)^(m + i*p), x], x] /; FreeQ[{a, b, i, m, p}, x] && IntegerQ[i] & 
&  !IntegerQ[p]
 

rule 53
Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int 
[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n, x], x] /; FreeQ[{a, b, c, d, n}, 
x] && IGtQ[m, 0] && ( !IntegerQ[n] || (EqQ[c, 0] && LeQ[7*m + 4*n + 4, 0]) 
|| LtQ[9*m + 5*(n + 1), 0] || GtQ[m + n + 2, 0])
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.10.23.4 Maple [A] (verified)

Time = 0.41 (sec) , antiderivative size = 136, normalized size of antiderivative = 1.04

method result size
gosper \(-\frac {\sqrt {c \,x^{2}}\, \left (b x +a \right )^{1+n} \left (-b^{3} n^{3} x^{3}-6 b^{3} n^{2} x^{3}+3 a \,b^{2} n^{2} x^{2}-11 b^{3} n \,x^{3}+9 a \,b^{2} n \,x^{2}-6 b^{3} x^{3}-6 a^{2} b n x +6 a \,b^{2} x^{2}-6 a^{2} b x +6 a^{3}\right )}{b^{4} x \left (n^{4}+10 n^{3}+35 n^{2}+50 n +24\right )}\) \(136\)
risch \(-\frac {\sqrt {c \,x^{2}}\, \left (-b^{4} n^{3} x^{4}-a \,b^{3} n^{3} x^{3}-6 b^{4} n^{2} x^{4}-3 a \,b^{3} n^{2} x^{3}-11 b^{4} n \,x^{4}+3 a^{2} b^{2} n^{2} x^{2}-2 x^{3} a n \,b^{3}-6 b^{4} x^{4}+3 a^{2} n \,x^{2} b^{2}-6 x \,a^{3} n b +6 a^{4}\right ) \left (b x +a \right )^{n}}{x \left (3+n \right ) \left (4+n \right ) \left (2+n \right ) \left (1+n \right ) b^{4}}\) \(156\)

input
int(x^2*(b*x+a)^n*(c*x^2)^(1/2),x,method=_RETURNVERBOSE)
 
output
-1/b^4/x*(c*x^2)^(1/2)*(b*x+a)^(1+n)/(n^4+10*n^3+35*n^2+50*n+24)*(-b^3*n^3 
*x^3-6*b^3*n^2*x^3+3*a*b^2*n^2*x^2-11*b^3*n*x^3+9*a*b^2*n*x^2-6*b^3*x^3-6* 
a^2*b*n*x+6*a*b^2*x^2-6*a^2*b*x+6*a^3)
 
3.10.23.5 Fricas [A] (verification not implemented)

Time = 0.24 (sec) , antiderivative size = 153, normalized size of antiderivative = 1.17 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=\frac {{\left (6 \, a^{3} b n x + {\left (b^{4} n^{3} + 6 \, b^{4} n^{2} + 11 \, b^{4} n + 6 \, b^{4}\right )} x^{4} - 6 \, a^{4} + {\left (a b^{3} n^{3} + 3 \, a b^{3} n^{2} + 2 \, a b^{3} n\right )} x^{3} - 3 \, {\left (a^{2} b^{2} n^{2} + a^{2} b^{2} n\right )} x^{2}\right )} \sqrt {c x^{2}} {\left (b x + a\right )}^{n}}{{\left (b^{4} n^{4} + 10 \, b^{4} n^{3} + 35 \, b^{4} n^{2} + 50 \, b^{4} n + 24 \, b^{4}\right )} x} \]

input
integrate(x^2*(b*x+a)^n*(c*x^2)^(1/2),x, algorithm="fricas")
 
output
(6*a^3*b*n*x + (b^4*n^3 + 6*b^4*n^2 + 11*b^4*n + 6*b^4)*x^4 - 6*a^4 + (a*b 
^3*n^3 + 3*a*b^3*n^2 + 2*a*b^3*n)*x^3 - 3*(a^2*b^2*n^2 + a^2*b^2*n)*x^2)*s 
qrt(c*x^2)*(b*x + a)^n/((b^4*n^4 + 10*b^4*n^3 + 35*b^4*n^2 + 50*b^4*n + 24 
*b^4)*x)
 
3.10.23.6 Sympy [F]

\[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=\begin {cases} \frac {a^{n} x^{3} \sqrt {c x^{2}}}{4} & \text {for}\: b = 0 \\\int \frac {x^{2} \sqrt {c x^{2}}}{\left (a + b x\right )^{4}}\, dx & \text {for}\: n = -4 \\\int \frac {x^{2} \sqrt {c x^{2}}}{\left (a + b x\right )^{3}}\, dx & \text {for}\: n = -3 \\\int \frac {x^{2} \sqrt {c x^{2}}}{\left (a + b x\right )^{2}}\, dx & \text {for}\: n = -2 \\\int \frac {x^{2} \sqrt {c x^{2}}}{a + b x}\, dx & \text {for}\: n = -1 \\- \frac {6 a^{4} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {6 a^{3} b n x \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} - \frac {3 a^{2} b^{2} n^{2} x^{2} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} - \frac {3 a^{2} b^{2} n x^{2} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {a b^{3} n^{3} x^{3} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {3 a b^{3} n^{2} x^{3} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {2 a b^{3} n x^{3} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {b^{4} n^{3} x^{4} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {6 b^{4} n^{2} x^{4} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {11 b^{4} n x^{4} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} + \frac {6 b^{4} x^{4} \sqrt {c x^{2}} \left (a + b x\right )^{n}}{b^{4} n^{4} x + 10 b^{4} n^{3} x + 35 b^{4} n^{2} x + 50 b^{4} n x + 24 b^{4} x} & \text {otherwise} \end {cases} \]

input
integrate(x**2*(b*x+a)**n*(c*x**2)**(1/2),x)
 
output
Piecewise((a**n*x**3*sqrt(c*x**2)/4, Eq(b, 0)), (Integral(x**2*sqrt(c*x**2 
)/(a + b*x)**4, x), Eq(n, -4)), (Integral(x**2*sqrt(c*x**2)/(a + b*x)**3, 
x), Eq(n, -3)), (Integral(x**2*sqrt(c*x**2)/(a + b*x)**2, x), Eq(n, -2)), 
(Integral(x**2*sqrt(c*x**2)/(a + b*x), x), Eq(n, -1)), (-6*a**4*sqrt(c*x** 
2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b**4*n**2*x + 50*b**4*n 
*x + 24*b**4*x) + 6*a**3*b*n*x*sqrt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10 
*b**4*n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24*b**4*x) - 3*a**2*b**2*n** 
2*x**2*sqrt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b**4*n 
**2*x + 50*b**4*n*x + 24*b**4*x) - 3*a**2*b**2*n*x**2*sqrt(c*x**2)*(a + b* 
x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24*b* 
*4*x) + a*b**3*n**3*x**3*sqrt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4* 
n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24*b**4*x) + 3*a*b**3*n**2*x**3*sq 
rt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b**4*n**2*x + 5 
0*b**4*n*x + 24*b**4*x) + 2*a*b**3*n*x**3*sqrt(c*x**2)*(a + b*x)**n/(b**4* 
n**4*x + 10*b**4*n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24*b**4*x) + b**4 
*n**3*x**4*sqrt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b* 
*4*n**2*x + 50*b**4*n*x + 24*b**4*x) + 6*b**4*n**2*x**4*sqrt(c*x**2)*(a + 
b*x)**n/(b**4*n**4*x + 10*b**4*n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24* 
b**4*x) + 11*b**4*n*x**4*sqrt(c*x**2)*(a + b*x)**n/(b**4*n**4*x + 10*b**4* 
n**3*x + 35*b**4*n**2*x + 50*b**4*n*x + 24*b**4*x) + 6*b**4*x**4*sqrt(c...
 
3.10.23.7 Maxima [A] (verification not implemented)

Time = 0.22 (sec) , antiderivative size = 116, normalized size of antiderivative = 0.89 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=\frac {{\left ({\left (n^{3} + 6 \, n^{2} + 11 \, n + 6\right )} b^{4} \sqrt {c} x^{4} + {\left (n^{3} + 3 \, n^{2} + 2 \, n\right )} a b^{3} \sqrt {c} x^{3} - 3 \, {\left (n^{2} + n\right )} a^{2} b^{2} \sqrt {c} x^{2} + 6 \, a^{3} b \sqrt {c} n x - 6 \, a^{4} \sqrt {c}\right )} {\left (b x + a\right )}^{n}}{{\left (n^{4} + 10 \, n^{3} + 35 \, n^{2} + 50 \, n + 24\right )} b^{4}} \]

input
integrate(x^2*(b*x+a)^n*(c*x^2)^(1/2),x, algorithm="maxima")
 
output
((n^3 + 6*n^2 + 11*n + 6)*b^4*sqrt(c)*x^4 + (n^3 + 3*n^2 + 2*n)*a*b^3*sqrt 
(c)*x^3 - 3*(n^2 + n)*a^2*b^2*sqrt(c)*x^2 + 6*a^3*b*sqrt(c)*n*x - 6*a^4*sq 
rt(c))*(b*x + a)^n/((n^4 + 10*n^3 + 35*n^2 + 50*n + 24)*b^4)
 
3.10.23.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 300 vs. \(2 (123) = 246\).

Time = 0.29 (sec) , antiderivative size = 300, normalized size of antiderivative = 2.29 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx={\left (\frac {6 \, a^{4} a^{n} \mathrm {sgn}\left (x\right )}{b^{4} n^{4} + 10 \, b^{4} n^{3} + 35 \, b^{4} n^{2} + 50 \, b^{4} n + 24 \, b^{4}} + \frac {{\left (b x + a\right )}^{n} b^{4} n^{3} x^{4} \mathrm {sgn}\left (x\right ) + {\left (b x + a\right )}^{n} a b^{3} n^{3} x^{3} \mathrm {sgn}\left (x\right ) + 6 \, {\left (b x + a\right )}^{n} b^{4} n^{2} x^{4} \mathrm {sgn}\left (x\right ) + 3 \, {\left (b x + a\right )}^{n} a b^{3} n^{2} x^{3} \mathrm {sgn}\left (x\right ) + 11 \, {\left (b x + a\right )}^{n} b^{4} n x^{4} \mathrm {sgn}\left (x\right ) - 3 \, {\left (b x + a\right )}^{n} a^{2} b^{2} n^{2} x^{2} \mathrm {sgn}\left (x\right ) + 2 \, {\left (b x + a\right )}^{n} a b^{3} n x^{3} \mathrm {sgn}\left (x\right ) + 6 \, {\left (b x + a\right )}^{n} b^{4} x^{4} \mathrm {sgn}\left (x\right ) - 3 \, {\left (b x + a\right )}^{n} a^{2} b^{2} n x^{2} \mathrm {sgn}\left (x\right ) + 6 \, {\left (b x + a\right )}^{n} a^{3} b n x \mathrm {sgn}\left (x\right ) - 6 \, {\left (b x + a\right )}^{n} a^{4} \mathrm {sgn}\left (x\right )}{b^{4} n^{4} + 10 \, b^{4} n^{3} + 35 \, b^{4} n^{2} + 50 \, b^{4} n + 24 \, b^{4}}\right )} \sqrt {c} \]

input
integrate(x^2*(b*x+a)^n*(c*x^2)^(1/2),x, algorithm="giac")
 
output
(6*a^4*a^n*sgn(x)/(b^4*n^4 + 10*b^4*n^3 + 35*b^4*n^2 + 50*b^4*n + 24*b^4) 
+ ((b*x + a)^n*b^4*n^3*x^4*sgn(x) + (b*x + a)^n*a*b^3*n^3*x^3*sgn(x) + 6*( 
b*x + a)^n*b^4*n^2*x^4*sgn(x) + 3*(b*x + a)^n*a*b^3*n^2*x^3*sgn(x) + 11*(b 
*x + a)^n*b^4*n*x^4*sgn(x) - 3*(b*x + a)^n*a^2*b^2*n^2*x^2*sgn(x) + 2*(b*x 
 + a)^n*a*b^3*n*x^3*sgn(x) + 6*(b*x + a)^n*b^4*x^4*sgn(x) - 3*(b*x + a)^n* 
a^2*b^2*n*x^2*sgn(x) + 6*(b*x + a)^n*a^3*b*n*x*sgn(x) - 6*(b*x + a)^n*a^4* 
sgn(x))/(b^4*n^4 + 10*b^4*n^3 + 35*b^4*n^2 + 50*b^4*n + 24*b^4))*sqrt(c)
 
3.10.23.9 Mupad [B] (verification not implemented)

Time = 0.37 (sec) , antiderivative size = 214, normalized size of antiderivative = 1.63 \[ \int x^2 \sqrt {c x^2} (a+b x)^n \, dx=\frac {{\left (a+b\,x\right )}^n\,\left (\frac {x^4\,\sqrt {c\,x^2}\,\left (n^3+6\,n^2+11\,n+6\right )}{n^4+10\,n^3+35\,n^2+50\,n+24}-\frac {6\,a^4\,\sqrt {c\,x^2}}{b^4\,\left (n^4+10\,n^3+35\,n^2+50\,n+24\right )}+\frac {6\,a^3\,n\,x\,\sqrt {c\,x^2}}{b^3\,\left (n^4+10\,n^3+35\,n^2+50\,n+24\right )}+\frac {a\,n\,x^3\,\sqrt {c\,x^2}\,\left (n^2+3\,n+2\right )}{b\,\left (n^4+10\,n^3+35\,n^2+50\,n+24\right )}-\frac {3\,a^2\,n\,x^2\,\sqrt {c\,x^2}\,\left (n+1\right )}{b^2\,\left (n^4+10\,n^3+35\,n^2+50\,n+24\right )}\right )}{x} \]

input
int(x^2*(c*x^2)^(1/2)*(a + b*x)^n,x)
 
output
((a + b*x)^n*((x^4*(c*x^2)^(1/2)*(11*n + 6*n^2 + n^3 + 6))/(50*n + 35*n^2 
+ 10*n^3 + n^4 + 24) - (6*a^4*(c*x^2)^(1/2))/(b^4*(50*n + 35*n^2 + 10*n^3 
+ n^4 + 24)) + (6*a^3*n*x*(c*x^2)^(1/2))/(b^3*(50*n + 35*n^2 + 10*n^3 + n^ 
4 + 24)) + (a*n*x^3*(c*x^2)^(1/2)*(3*n + n^2 + 2))/(b*(50*n + 35*n^2 + 10* 
n^3 + n^4 + 24)) - (3*a^2*n*x^2*(c*x^2)^(1/2)*(n + 1))/(b^2*(50*n + 35*n^2 
 + 10*n^3 + n^4 + 24))))/x